Fermi velocity engineering in graphene by substrate modification 
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The Fermi velocity, vf, is one of the key concepts in the study of a material, as it bears informa- 
tion on a variety of fundamental properties. Upon increasing demand on the device applications, 
graphene is viewed as a prototypical system for engineering vf- Indeed, several efforts have suc- 
ceeded in modifying vp by varying charge carrier concentration, n. Here we present a powerful 
but simple new way to engineer vf while holding n constant. We find that when the environment 
embedding graphene is modified, the vf of graphene is (i) inversely proportional to its dielectric 
constant, reaching Wp ~2.5xf0 6 m/s, the highest value for graphene on any substrate studied so 
far and (ii) clearly distinguished from an ordinary Fermi liquid. The method demonstrated here 
provides a new route toward Fermi velocity engineering in a variety of two-dimensional electron 
systems including topological insulators. 



I. INTRODUCTION 



while keeping n constant. 



Due to its lattice structure and position of the Fermi 
energy, the low energy electronic excitations of graphene 
are described by an effective field theory that is Lorentz 
invariant Unlike Galilean invariant theories such as 
Fermi liquids whose main relevant parameter is the ef- 
fective mass, Lorentz invariant theories are characterized 
by an effective velocity. Because of this, an increase of 
electron-electron interactions induces an increase of the 
Fermi velocity, vp, in contrast to Fermi liquids, where the 
opposite trend is true @. In the case of graphene, when 
electron-electron interactions are weak Q , i>f is expected 
to be as low as 0.85xl0 6 m/s, whereas, for the case of 
strong interactions [5|, vp is expected to be as high as 
1.73xl0 6 m/s. 

Recently, Fermi velocities as high as ~3x 10 6 m/s 
have been achieved in suspended graphene through a 
change of the carrier concentration n [fHH- However, 
because this dependence is logarithmic, n needs to be 
changed by two orders of magnitude in order to change 
the velocity by a factor of 3. This implies that it is un- 
practical to use n as a way to engineer vp, let alone the 
fact that one should first realize suspended graphene in 
the device @. Several other routes have also been pro- 
posed to engineer vp in graphene via the electron-electron 
interaction, including modifications of: a) curvature of 
the graphene sheet Ejjhb) periodic potentials c) 
dielectric screening |12h14| . While the former two also 
substantially modify the starting material, the latter sim- 
ply modifies the effective dielectric constant, e, making it 
more appealing for device applications 15j. Despite this 
advantage, no systematic study of how to engineer vp by 
changing e exists to date. Here we provide a new venue to 
control the Fermi velocity of graphene using dielectrics, 
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II. RESULTS 

We perform such a study using three single-layer 
graphene samples, which were prepared by chemical va- 
por deposition (CVD) on Cu, followed by an in situ 
dewetting of Cu on quartz (single crystal SiO ?) fig or 
a transfer onto hexagonal boron nitride (BN) [17j . and 
by epitaxial growth on Lff-SiC(OOOl) p]|. Figures 1A 
and IB show angle- resolved photoemission spectroscopy 
(ARPES) intensity maps measured near the Brillouin 
zone (BZ) corner K along the T-K direction for the two 
CVD grown samples, which constitute the first report 
on Dirac quasiparticle mapping from these samples. Fol- 
lowing the maximum intensity, one can clearly observe 
almost linear energy spectra, characteristic of Dirac elec- 
trons 19]. The momentum distribution curves (MDC), 
intensity spectra taken at constant energy as a function 
of momentum, are shown in Fig. 1C. In addition to be- 
ing proportional to the imaginary part of the electron 
self-energy, the MDC spectral width provides informa- 
tion on the sample quality. A clear increase of the width 
is observed by changing the substrate from SiC(OOOl) 
via BN to quartz, a trend that is in overall agreement 
with the theoretical expectation that the electron self- 
energy should vary with the inverse square of the dielec- 
tric screening [2(| , as later discussed. The quartz sample 
here used constitutes a substantial improvement over a 
previous experiment [21( on a similar substrate (compare 
0.19 A" 1 (red line) versus ~0.7 A 1 (gray-dashed line)). 
The much improved data quality allows for a detailed 
self-energy analysis and consequent extraction of impor- 
tant parameters such as vp. 

To understand how the dielectric substrate affects the 
electronic properties, in Fig. 2, we show the energy vs. 
momentum dispersions for graphene on three different 
substrates, SiC(OOOl), BN, and quartz, obtained by fit- 
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FIG. 1. ARPES intensity maps of graphene on quartz and 
BN. (A-B) Normalized and raw ARPES intensity maps of 
graphene/quartz (panel (^4)) and graphene/BN (panel (B)), 
respectively. The red and dark-yellow lines are the dis- 
persions, obtained by fitting momentum distribution curves 
(MDCs). (O MDCs at E F for graphene on SiC(000T) (blue 
line), BN (dark-yellow line), quartz (red line), and Si02 [H[ 
(gray-dashed line). 



FIG. 2. Experimental and theoretical energy spectra for dif- 
ferent dielectric constants. (A) Experimental dispersions for 
graphene on SiC(0001) (blue line), BN (dark- yellow line), and 
quartz (red line). The gray-dashed lines are guides to the 
eyes. The insets are cartoons for the electron band struc- 
ture of graphene with weak (left) and strong (right) electron- 
electron interactions. The data are shifted along the x-axis. 
(B) The direct comparison of experimental dispersions with 
theories: e = oo (magenta line) [|[ and e = 1 (cyan line) [T3 |. 



ting the MDC spectra. The observed dispersions exhibit 
two distinctive features. First, the measured dispersions 
deviate from linearity with an increased slope around 
~ —0.5 eV for all the samples (compare experimental 
data to dashed gray lines in Fig. 2 A). As the substrate is 
changed from SiC(0001) via BN to quartz, correspond- 
ing to a decrease of the dielectric screening, the departure 
from linearity at high energy becomes more pronounced. 
Second, the direct comparison between experimental dis- 
persions and ab initio calculations for the two extreme 
cases e = 1 [1| (suspended graphene) and e = oo [|| shows 
another substrate-dependence (Fig. 2B). Upon changing 
the substrate, the slope increases approaching the dis- 
persion for e = 1. The deviation from linearity and the 
enhancement of the slope result in a reshape of the typi- 
cal conical dispersion, in a similar fashion as reported for 
other charge-neutral graphene samples @, [l2j (see car- 
toons in the inset of Fig. 2 A: from left to right). We 
note that the largest upturn for graphene/quartz cannot 
be explained by: a) resolution, which typically results in 
the deflection of MDC peaks near Ep to lower momen- 
tum, and would involve a much smaller effect by an order 
of magnitude (<a few tens meV) [22| ; b) the presence of 
other bands with a different azimuthal orientation, which 
would cause instead an abrupt increase and a significant 
asymmetry of the MDC width at the upturn energy. 



III. DISCUSSION 

To quantify the effect of dielectric substrates on the 
electron-electron interactions and vp , we adopt the stan- 
dard self-energy analysis to extract self-consistently the 
strength of the electron-electron interactions and e [l], [TH, 
US Hil- Figure 3^4 shows the difference between measured 
dispersions, E(k) (from Fig. 2 A), and the theoretical dis- 



persion for e = oo, i?LDA(k) (shown in Fig. 2B). As- 
suming that electron-electron interactions are effectively 
screened for e = oo, the E — E'lda curve can be con- 
sidered a good measurement of the difference between 
the self-energy and its value at Ep. To fit these curves, 
we use the marginal Fermi liquid self-energy function 
as previously reported [12l . |23[ with an analytic form of 
^ (k-k F ) lnj^ (dotted lines in Fig. 3 A). Here, a is 
a dimensionless fine-structure constant (or the strength 
of electron-electron interactions) defined as iT . e ehv [23j j . 
vq the Fermi velocity for e = oo, 0.85 xlO 6 m/s ||, k c 
the momentum cut-off, 1.7 A^ 1 , and kp the Fermi wave 
number. An overall good agreement with the experi- 
mental data is observed allowing to extract important 
parameters such as e and a for graphene on each sub- 
strate. For graphene on SiC(0001) and BN, we obtain 
e = 7.26±0.02 (a = 0.35) and e = 4.22±0.01 (a = 0.61), 
respectively. The extracted value for graphene on BN 
is in agreement with the standard approximation e = 
(evacuum+e S ubstratc)/2 = 4.02 and 3.05, where e vacnum = 1 
and Csubstrato = 7.04 (for out-of-plane polarization) and 
5.09 (for in-plane polarization) in the low frequency limit 
(static dielectric constant) for hexagonal-BN [25(. Sim- 
ilarly, the obtained value f or g raphene on SiC(0001) is 
close to a previous report [12} ■ The apparent discrep- 
ancy with the latter (compare e = 7.26±0.02 in this 
work with 6.4±0.1 in reference 12]) is due to the dif- 
ferent choice of reference band (or so-called bare band). 
Specifically, in this work, -Elda is used as the bare band, 
whereas, in reference [12j, the bare band is approximated 
by a straight line. Finally, for graphene/quartz, we ob- 
tain e = 1.80±0.02 (a — 1.43), which is smaller than 
the expected value of e = 2.45 [26|], instead closer to the 
experimentally extracted value for suspended graphene 
(-2.2) @. This observation, together with the simi- 
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lar energy-momentum dispersion relation at high bind- 
ing energy to the theoretical one for suspended graphene 
(Fig. 25), point to a very weak effect of the substrate. 
This is likely a consequence of the different sample prepa- 
ration method adopted here (see Materials and Methods 
section). 

In Fig. 35, we show the measured vp as a function 
of the extracted e (see also Table 1). Results from a 
suspended sample [6| and another graphene/Si02 sam- 
ple are also plotted for comparison. Upon de- 
creasing e from oo to 7.26 and 4.22, vf is enhanced 
from its LDA limit of 0.85 xlO 6 m/s (cyan triangle in 
Fig. 35) to I.f5±0.02xl0 6 m/s (blue circle in Fig. 35) 
and 1.49±0.08xf0 6 m/s (dark- yellow circle in Fig. 35), 
by 35 % and 75 %, respectively. Surprisingly, when e 
is further decreased to 1.80, a dramatic enhancement of 
v F up to 2.49±0.30xl0 6 m/s (red circle in Fig. 35) is 
observed. Such enhancement corresponds to a 190 % in- 
crease from its bare value and represents the hig hest value 
reported for graphene on any substrate |27H29j | . Interest- 
ingly, this velocity is comparable to the value measured 
for suspended graphene (green square in Fig. 35) [6j. 
Clearly, a 1/e dependence of vp is observed (dashed line 
in Fig. 35) in agreement with the theoretical predic- 
tion [a, [HI- Our result constitutes the first observation of 
a power law dependence of the Fermi velocity on the di- 
electric constant at fixed n. This power law dependence 
allows one to achieve, by a smart choice of dielectric, a 
high value of vf that cannot be attained otherwise by 
changing n [6J. 

We note that CVD graphene on quartz (red circle 
in Fig. 35) exhibits higher vp than exfoliated graphene 
on amorphous SiC>2 (gray square in Fig. 35) with the 
same stoichiometry as quartz. This is a consequence of 
different sample preparation process and is due to the 
larger presence of impurities in the exfoliated sample, 
as suggested by the extremely broad spectra (see gray 
dashed line in Fig. \C). Therefore, although, in theory, 
one should expect smaller vp due to screened electron- 
electron interactions from impurity [13], one should be 
cautious in extracting meaningful parameters from these 
data. We also note that ab initio GW calculations [5] 
(magenta triangle in Fig. 35) underestimate vf of sus- 
pended graphene. This may be due to the finite fc-point 
sampling inherent in such calculations, or it could also 
be an indication of the need to add higher-order terms in 
the self-energy calculation by the GW-approximation. 



TABLE I. Fermi velocity (vf), dielectric constant (e), and fine 
structure constant (a) of graphene on each substrate 



Substrate 


vfxIO 6 m/s 


c 


a 


Metals (LDA) 


0.85 


oo 




SiC(OOOl) 


1.15±0.02 


7.26±0.02 


0.35 


fe-BN 


1.49±0.08 


4.22±0.01 


0.61 


Quartz 


2.49±0.30 


1.80±0.02 


1.43 



In Fig. 3C, we plot the ratio between vf and vq, the ex- 
pected Fermi velocity in the fully screened case (e = oo), 
as a function of a. As the strength of electron-electron 
interactions is increased, vf is also enhanced. This is in 
striking difference with the standard Fermi liquid picture, 
where vp is expected to decrease with increasing a [30j ] . 
On the other hand, the observed behavior is consistent 
with previous theoretical studies for graphene in the case 
of specific electron-electron interactions [3(| HH (dashed 
line in Fig. 3(7) exhibiting the characteristic self-energy 
spectrum analogous to a marginal Fermi liquid [l|. As 
a result, the departure from the Fermi liquid picture be- 
comes more important with increasing electron-electron 
interactions or decreasing dielectric screening (see the re- 
lation between a and e in the inset of Fig. 3C). Addition- 
ally, the observation of a values close to 1 (neither a <C 1 
nor a>l) for graphene/quartz may indicate that a full 
theoretical treatment beyond the random-phase approx- 
imation [l| may be required to understand this sample 
and/or suspended graphene @. 

The very good agreement with theoretical predic- 
tions (23|, |3l| for both vp versus e (Fig. 35) and vp ver- 
sus a (Fig. 3(7) confirms that the dielectric constants 
obtained by the self-energy analysis are self-consistent. 
Finally the experimentally determined e can largely 
account for the relatively broad MDCs observed for 
graphene on quartz (Fig. 1(7), as compared to graphene 
on BN and SiC(0001). For e values of 1.80, 4.22, and 
7.26, for graphene on quartz, BN, and SiC(0001) respec- 
tively, the MDC widths, expected to vary with the inverse 
square of the dielectric screening ;20] , should be roughly 
16 and 5 times broader for graphene on quartz and BN 
than graphene on SiC(0001), in line with the experimen- 
tal observation (see, for example, Fig. 1(7). We stress 
that, contrary to a Fermi liquid system, the broader MDC 
spectra observed for graphene/quartz do not necessarily 
imply decreased transport properties. On the contrary, 
the enhanced a, the primary cause of the broad spec- 
tra, give rise to an enhancement of Fermi velocity, which 
is ultimately one of the most important parameters for 
device applications. 

In conclusion, we have unveiled the crucial role of di- 
electric screening in graphene to control both Fermi ve- 
locity and electron-electron interactions. Additionally, 
we have shown that graphene, in its charge neutral state, 
departs from a standard Fermi liquid not only in its log- 
arithmic energy spectrum as previously discussed [12j . 
but also in the way that vp is modulated by the strength 
of electron-electron interactions. This dependence pro- 
vides an alternative way to engineer Fermi velocity for 
graphene on a substrate by modifying the dielectric sub- 
strate. This approach can also be applied to charge- 
doped graphene and other two-dimensional electron sys- 
tems such as topological insulators [32| that can be grown 
or transferred to dielectric substrates. 
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FIG. 3. Fermi velocity and the strength of electron-electron interactions. (A) E — £"lda dispersions for graphene on SiC(OOOl) 
(blue line), BN (dark- yellow line), and quartz (red line). (B) Fermi velocities as a function of e. The dashed line is a theoretical 
curve for vf, which is inversely proportional to e [HHH]. Filled symbols correspond to experimental results, while empty symbols 
to theoretical values, e = 2.45 for G/Si02 [IB] is obtained from the standard approximation, e = (e V acuum + £substrate)/2 (see 
text). (C) The ratio of vf, the renormalized Fermi velocity due to electron-electron interactions, to vo = 0.85xl0 6 m/s, the 
bare Fermi velocity in the LDA limit where e = oo Q, as a function of a. The dashed line is the fit given by vf/vo = 
1 — 3.28a{l + (l/4)ln[(l + 4a)/4a] — 1.45} [31| for charge neutral graphene. The inset is the relation between a and e, where 
the dashed line is a a = in ^ h = 2.57/e curve [2^ ]. 
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IV. MATERIALS AND METHODS 

Graphene samples were prepared in three different 
ways: epitaxial growth on the surface of a 4ii-SiC(0001) 
substrate; chemical vapor deposition (CVD) growth on a 
Cu film followed by a transfer onto the surface of boron 
nitride [171 ]: and CVD growth followed by in situ dewet- 
ting of Cu layer in between graphene and a single crys- 
tal Si02 (namely quartz which is different from amor- 
phous Si02 on an Si substrate, the widely used substrate 
for exfoliated graphene 27J) substrate The later 

procedure is clearly different from the standard method 
of exfoliating graphite followed by deposition onto the 
amorphous S1O2 layer [2l|. This results in a reduced ef- 



fect of the substrate that is suggested by the enhanced 
height variation with respect to the substrate compared 
to the s amp le prepared by the exfoliation and deposi- 
tion [ll| [33 ■ The resulting graphene is more decoupled 
from the substrate as supported by several features such 
as Fermi velocity, dielectric constant, and the electron 
band at higher energies closer to suspended sample. 

In order to remove any residue including Cu and 
PMMA, a precursor to grow CVD graphene and a poly- 
mer to transfer graphene, respectively, we heated the 
sample to 1000 °C in ultra- high vacuum. The removal 
of Cu is confirmed by: (a) optical microscopy showing a 
cleaner image without residual Cu once the sample has 
been heated; (b) absence of related Cu features in the 
ARPES spectra such as 3d electrons at 3.0 eV and 3.5 eV 
below Fermi energy, and 4s free-electron-like state with 
a band minimum at 0.25 eV below Fermi energy pjij . 

High-resolution ARPES experiments have been per- 
formed at beamline 10.0.1.1 of the Advanced Light 
Source at Lawrence Berkeley National Laboratory using 
50 eV photons at 15 K. Energy and angular (momentum) 
resolutions were set to be 22 meV and 0.2 (^0.01 A" 1 ), 
respectively. 
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